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Abstract 

It is known that all spatially homogeneous solutions of the vacuum 
Einstein equations in four dimensions which exist for an infinite proper 
time towards the future are future geodesically complete. This paper 
investigates whether the analogous statement holds in higher dimensions. 
A positive answer to this question is obtained for a large class of models 
which can be studied with the help of Kaluza-Klein reduction to solutions 
of the Einstein-scalar field equations in four dimensions. The proof of 
this result makes use of a criterion for geodesic completeness which is 
applicable to more general spatially homogeneous models. 

1 Introduction 

The solutions of the Einstein equations most commonly applied in cosmology are 
spatially homogeneous and include some that recollapse and some that expand 
for ever. In the former case the proper length of the timelike curves orthogonal 
to the hypersurfaces of homogeneity (the worldlines of comoving observers) is 
finite. Along a curve of this type the expansion of this congruence of curves is 
first positive, then zero and finally negative. This is what is meant by recollapse. 
In the latter case the expansion is always positive and the length of the curves 
orthogonal to the hypersurfaces of homogeneity is infinite towards the future. 
This is what is meant by saying that the solution expands for ever. 

On physical grounds it is reasonable to expect that a forever expanding 
solution will be future geodesically complete, i.e. there should be no singularities 
in the future. A theorem which confirms this expectation is proved in [17j 
(Theorem 2.1 of that paper). In order to explain this result let t be a Gaussian 
time coordinate based on one of the hypersurfaces of homogeneity and let trfc 
be the mean curvature of these hypersurfaces. Then trfc is a function of t. The 
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expansion of the congruence of timelike geodesies mentioned above is — |trfe. 
Define an expanding phase of the solution to be a time interval (ti, ia) on which 
trfc < 0. Here t\ and may be finite or infinite. The theorem says that, 
under certain assumptions, if t 2 < 00 the solution can be extended to a time 
interval (£1,^3) with £3 > t^. Moreover, if ti = 00 then the spacetime is future 
geodesically complete. The theorem is formulated for spatially homogeneous 
solutions of the Einstein equations coupled to a matter model satisfying three 
assumptions. Two of these are inequalities on the energy-momentum tensor, 
the dominant energy and non-negative sum pressures conditions. The third is a 
continuation property for solutions of the Einstein- matter equations. The only 
thing about these conditions which is important in what follows is that they are 
all satisfied in the case of the Einstein vacuum equations. Hence the theorem 
applies in the vacuum case. 

In a spatially homogeneous spacetime, by definition the isometry group acts 
transitively on the hypersurfaces of homogeneity. In four dimensions there is a 
well-known classification of the possible isometry groups. Either the group can 
be taken to have dimension three or a higher-dimensional group is necessary. In 
the second of these cases there is only one possibility, which is called Kantowski- 
Sachs symmetry. In the first case any connected three-dimensional Lie group 
can occur. For the purposes of what follows it may be assumed without loss of 
generality that the spacetime is simply connected. This is because passing to 
the universal cover of a given spacetime has no effect on the dynamics. With 
this assumption the hypersurfaces of homogeneity can be identified with the 
group itself and the induced metric on one of these hypersurfaces with a left- 
invariant metric on that group. The classification of connected and simply 
connected Lie groups reduces to that of the corresponding Lie algebras. The 
three-dimensional Lie algebras were classified by Bianchi into types I to IX. It 
turns out that a maximally extended vacuum spacetime of Bianchi type I- VIII 
expanding at some time is forever expanding while one which is of type IX or has 
Kantowski-Sachs symmetry recollapses [14]. The cases where recollapse takes 
place are precisely those where there is a spatial metric in the class whose scalar 
curvature R is positive. That this is a necessary condition for recollapse is not 
hard to see. The Hamiltonian constraint on a hypersurface of homogeneity is 
(in the vacuum case) 

R - k ab k ab + (trfc) 2 = (1) 

where k a b is the second fundamental form. Recollapse means that there is a 
maximal hypersurface, i.e. one where trfc = 0. There R is clearly non- negative. 
If it is zero then k a b = 0. Moreover it can be shown that if R = the induced 
metric of that hypersurface is flat and that the solution arising from those data is 
a flat spacetime where trfc vanishes at all times. The latter case cannot occur in 
the context of spatially homogeneous cosmological models which are expanding 
at some time. 

The aim of this paper is to investigate to what extent the results just stated 
for four-dimensional spatially homogeneous vacuum spacetimes extend to higher 
dimensions n + 1. In particular, sufficient conditions for future geodesic com- 
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pleteness are derived. Only the case where the isometry group can be taken 
to be of dimension n is considered. This is the analogue in higher dimensions 
of Bianchi models. The analogue of Kantowski-Sachs models is not considered. 
The same argument as in four dimensions applies to show that a necessary con- 
dition for recollapse is the existence of a left- invariant metric of positive scalar 
curvature on the group of interest. This motivates us to restrict attention to 
Lie groups admitting no left-invariant metric of positive scalar curvature. In 
this case it has been shown that vacuum spacetimes with this type of symmetry 
exist for infinite Gaussian time in the future ([T8], Theorem 5.3). On the other 
hand it turns out that the arguments used to prove geodesic completeness in 
the four-dimensional case do not easily generalize to higher dimensions. In what 
follows partial results on this question are obtained. 

For a general dimension there is no known analogue of the Bianchi classifica- 
tion of Lie algebras. In some of the lower dimensions there are results. A review 
of what is known in the case of four-dimensional Lie algebras can be found in [8] . 
The strategy in the following is to identify some classes of Lie groups for which 
geodesic completeness can be proved. The main case considered is that where 
there is an isometry group of the form G3 x (R xi Z2)™ -3 for a three-dimensional 
group G3. The spatial manifold is given by the identity component of this group, 
G 3 x K n ~ 3 . The strategy for handling these groups is to use Kaluza-Klein reduc- 
tion. The presence of the discrete symmetry defined by the group Z^~ 3 ensures 
that the result is the Einstein equations coupled to linear scalar fields in four 
dimensions. Without it a more complicated reduced system would be obtained. 
In the case n = 4, for instance, it would be the Einstein-Maxwell-dilaton sys- 
tem. Section 2 contains the basic notation and equations, including those for 
the Kaluza-Klein reduction. In the third section results on the Einstein-scalar 
field system in four dimensions are proved. These are used in treating higher- 
dimensional vacuum spacetimes in later sections and are also of interest in their 
own right. A general criterion for geodesic completeness is stated and proved 
in Section 4. It is then combined with the results on the scalar field to obtain 
the main theorem, Theorem \5\ Section 5 discusses geodesic completeness in 
some examples with other types of symmetry and illustrates how the results of 
Section 4 can be applied. Conclusions and an outlook are contained in Section 
6. 

This paper is based in part on the diploma thesis of the first author [7]. 

2 The basic equations 

Let G be a connected and simply connected Lie group and 9 a a basis of the 
space of left-invariant one-forms on G. Consider the metric 

- dt 2 + g ab {t)9 a ® 6 b (2) 

on the manifold [to, 00) x G for some constant to- Let k a b be the components in 
the given basis of the second fundamental form of the level hypersurfaces of t. 
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By the definition of the second fundamental form 

d t g a b = -2fc a 6- (3) 

It will now be supposed that this metric satisfies the Einstein vacuum equations 
and these will be decomposed using the foliation by the level hypersurfaces of 
t. The Einstein constraint equations are 

R - k ab k ah + (trfc) 2 = 0, (4) 

V a k a b = 0. (5) 

The evolution equations are given by 

d t k a b = (trfc)fc% + R\. (6) 

A real number A is said to be an eigenvalue of k ab with respect to g ab if there 
is a non-zero vector v a which satisfies the equation k ab v b — Xg ab v b . In n space 
dimensions there are n eigenvalues of k ab with respect to g ab , counting multi- 
plicity. Call them A a . The mean curvature trfc = ^ a A a never vanishes in an 
expanding phase. Define the generalized Kasner exponents (GKE) by p a = 
When the Lie group G is Abelian the p a are independent of time and the solution 
can be written explicitly as 

n 

-dt 2 +J2t 2pa {dx a ) 2 (7) 

a=l 

where Y^a=iPa — 1 and J2a=iPa = 1- I R tnc casc n = 3 this is the Kasner 
solution. 

These equations will now be specialized to the case that the Lie group is of 
the form G3 x R n_3 for a three-dimensional group G3 and the left-invariant basis 
9 a is adapted to the product decomposition. Call the basis vectors tangent to the 
first factor 9 1 , i = 1,2,3. The basis vectors corresponding to the second factor 
commute and so can be chosen to be of the form dy 1 , I = 1, . . . , n — 3. Suppose 
further that the action of the group extends to an action of G3 
which preserves the hypersurfaces of constant t and the product decomposition 
of the initial hypersurface. Then the generators of Z^ -3 act by commuting 
reflections of ]R™~ 3 . Without loss of generality they may be assumed to be 
given by reflections in the coordinates y 1 . The invariance of the metric under 



these transformations implies that it can be written in the form 

- dt 2 + 9lJ {t)6 l ®6 j e 20 ' ( 'W) 2 (8) 

i 

for scalar functions (p 1 . Let <f> — ^ 7 (j) 1 and write gij = e^^gij. Let t be 
a Gaussian coordinate for the conformally rescaled metric. Then = 
Rewriting the Einstein equations in terms of and (j) 1 gives the equations 

-kij~k ij + (trfc) 2 = p, (9) 
di¥j = (trfc)^- + - P 5) , (10) 
<9 t V - (trfc)<9 t ~/ = 0. (11) 
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Here 

p=\T,^) 2 - ( 12 ) 

These equations are identical to the Einstein equations in four dimensions cou- 
pled to n — 3 non-interacting massless scalar fields. For prescribed initial data 
it is possible to do a rotation in R™ -3 so as to set all d^y 1 except d^tj) 1 to 
zero. By the equations of motion for the scalar fields all cj) 1 except cj) 1 will be 
constant during the evolution. Then only cf) 1 makes a contribution to the energy- 
momentum tensor. The essential dynamics is that of a single scalar held. The 
four-dimensional metric and the scalar held dehne a five-dimensional metric. 
The full spacetime is obtained by taking the product of the five-dimensional 
metric with a fiat time-independent metric. 

The generalized Kasner exponents of the (n + l)-dimensional metric are 
related to the four-dimensional quantities by 

A a + \e-$d t <j> , 1Q . 

Pa = — _ ± , (13) 

trfc + 2 dtp 



2< 



-e 2 dt^ 1 



PI = r rrfrr - ( 14 ) 

trfc + 2 d t 4> 

In a spatially homogeneous situation a scalar held is equivalent to an untilted 
stiff fluid. Thus to understand the dynamics of the scalar field case known 
results for a stiff fluid can be applied. The matter variables of the two models 
are related by p — ^(d^c/)) 2 . Note that e~?d t = <3j. 

It was pointed out in the introduction that it is important whether a given 
Lie group G admits left-invariant metrics with positive scalar curvature. If G 
is a Lie group of any of the Bianchi types except IX it admits no metric of 
this kind. This is equivalent, for simply connected Lie groups, to the property 
that they are diffcomorphic to Euclidean space [2]. It follows that a Lie group 
with topology G x M™ -3 is also homeomorphic to Euclidean space and thus also 
admits no left-invariant metrics of positive scalar curvature. Hence groups of this 
type fail the necessary condition for recollapse mentioned in the introduction. 
This is a sign that they are good candidates for the symmetry groups of future 
geodesically complete metrics. 



3 Four-dimensional models with a scalar field 

This section contains results on the late-time behaviour of solutions of the 
Einstein-scalar field system in four dimensions, some of which are deduced from 
results on untilted stiff fluids taken from the literature. When information on 
the dynamics of the energy density is available the asymptotics of (f> can be 
obtained by integration. Much of the analysis is based on the Wainwright-Hsu 
formulation of the equations for Bianchi models [23], which will now be recalled. 
Since the only case of interest in the following is that of a stiff fluid the equations 
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are only written down for the case where the parameter 7 in the equation of 
state p = (7 — l)p of the fluid is equal to two. The system of evolution equations 
is 



N[ = (q-4Z + )Ni, (15) 

N' 2 = ( g + 2S+ + 2V3£_)7V 2 , (16) 

JV£ = (g + 2S+-2\/3E_)^3, (17) 

E' + = -(2-g)E + -3S+, (18) 

E'_ = -(2-g)E_ -3S1, (19) 

ft'=-2(2-g)ft (20) 



where 



? = 2(ft + E + +E 5 :), (21) 
S+ = ~[(JV 3 - iV 2 ) 2 - JVi(2JVi - 7V 2 _ jv 3 )], (22) 

5_ = -(Na-^XN^-Ni-Na). (23) 

The prime denotes -4p and t is related to the Gaussian time t by % — — |trfe. 
There is also the equation which expresses the Hamiltonian constraint in these 
variables which is 



ft + Y?+ + E 2 . + ^ + #2 + ^3 - 2(iV!^ 2 + 7V 2 7V 3 + N1N3)) = 1. (24) 



It will not be necessary to recall the definition of all the variables in the system 
here. What is important is that £+ and £_ are linear combinations of the 
generalized Kasner exponents which are linearly independent. This means that 
if £+ and £_ tend to finite limits as t — > 00 the same is true of pi, pi and 
p 3 . The density parameter ft is equal to , tr £ 2 . The mean curvature can be 
recovered from a solution of the Wainwright-Hsu system using the equation 

(trfc)' = -(l + <?)(trfc). (25) 

Theorem 1 Let (M, g) be a solution of the Einstein equations coupled to a 
scalar field <f> in four dimensions with symmetry of a Bianchi type other than 
IX or VI_i. If the solution is expanding at some time and the time interval of 
definition of the solution is maximal then the upper limit of the interval is infinity 
and the generalized Kasner exponents and the density parameter ft converge to 
limits as t — > 00. 



Proof Bianchi type I solutions are defined by the condition Ni = N 2 = N 3 = 0. 
It follows from the Hamiltonian constraint that q — 2 and that E+, S_ and 
ft are constant. That the conclusion of the theorem holds for Bianchi type II 
solutions, which are defined by the conditions N2 = iV 3 = and Ni ^ 0, was 
shown in Lemma 7.1 of [19]. Bianchi type VIq solutions are defined by the 
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conditions N± — 0, N2 > and 7V3 < 0. The evolution equation for E + reduces 
to 

E^ = -(2-g)(E+ + l). (26) 

It follows from the Hamiltonian constraint that q < 2. This in turn implies 
that £+ + 1 is strictly positive and that £ + is strictly decreasing. Applying 
the monotonicity principle to £+ shows that there can be no w-limit point with 
q < 2. Since the solution stays in a compact region it follows that q — > 2 as 
t — > 00 and that N\ and N2 tend to zero. Since £+ and il are monotone and 
bounded both of them tend to limits as r — > 00. By the Hamiltonian constraint 
the same is true of . Bianchi type VIIo solutions are defined by the conditions 
Ni — 0, N2 > and N3 > 0. Locally rotationally symmetric solutions of type 
VIIo are defined by the additional conditions N2 = N3 and £_ = 0. In this 
case S+, Si are zero and q = 2. Hence ft and £+ are constant. The proof 
for solutions of type VIIo which are not locally rotationally symmetric results 
from a small modification of the proof of Proposition 5 of [20], which treats 
the corresponding question in the vacuum case. First it can be shown by an 
argument similar to that used for type VIo case that and £+ are strictly 
monotone. This makes use of the equation 

{N 2 - N 3 )' = (q + 2£+)(V 2 - N 3 ) + 2%/3E_(V 2 + N 3 ) (27) 

to treat the points where E' + vanishes. Since fi and E + are also bounded they 
must converge as t — > 00. It is not known a priori that the whole solution is 
bounded and so it may not have an cj-limit point. Suppose that it does have such 
a point. Using the monotonicity principle it can be shown that if t„ is a sequence 
of times along which the solution converges then either + (A^2 — A^3) 2 )(r rl ) 
tends to zero or Ni + N% does so. It can then be concluded as in [20] . using 
the monotone function Z-i of that paper that >0asr— >oo. Note that 
Z_i is monotone in the case with matter (see [5D], Lemma 10.1). It remains to 
treat the case in which the solution has no uj limit point. Then N\ and N2 must 
tend to infinity and it can be assumed without loss of generality that the limit 
of SI + is less than one. In particular this means that the limit of £+ cannot 
be minus one. Hence, by (f2"6"]l the quantity 2 — q must be in L . It can, however, 
be shown by a slight modification of an argument of Ringstrom ([2D], proof of 
Proposition 5) that 2 — q is not in L . That argument uses two quantities called 
x and y and the necessary modification is to replace them by 



(28) 

%/3 ^2-^3 



(l-si-n)i 



2 (i-E5.-fi)i 



(29) 



Bianchi type VIII solutions are defined by the conditions iVi < 0, N2 > and 
N 3 < 0. The result of the theorem in this case follows from Theorem 3.1 and 
Corollary 3.1 of [12]. 
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For models of Bianchi class B it follows from Proposition 5.3 of [TT] that the 
limit set of any solution lies in a set The parameter k distinguishes the 
different Bianchi types and is equal to hr 1 in types VI?,, and VII^. For each 
value of fc, the set Ck consists of equilibrium solutions. It is an open curve 
Ck together with its two endpoints. The oj-limit set of a solution is in fact a 
single point of Ck- This can be seen as follows. A generic point of the curve 
is transversely hyperbolic and so if it belongs to the w-limit set no other point 
can. The remaining points form a discrete set and so by connectedness of the 
uj- limit set it cannot contain more than one of them. 

The fact that Bianchi type IX is excluded in the hypotheses of this theorem is 
essential. For the metric product of a four-dimensional solution of Bianchi type 
IX with the real line fails to exist globally towards the future. On the other 
hand there is no reason to expect that type VT_i has to be excluded. It appears 
nevertheless that up to now no proof is availabfe in that case. 

Some more precise information will now be given on the asymptotics of scalar 
field models. Let T> be the disk in the space with coordinates (£+, E_, Nx,N2, A3) 
defined by the conditions that the Ni are all zero and E+ + T?_ < 1. Let C be 
the line defined by N 1 = 0, A 2 = A 3 > 0, E+ = -1 and E_ = 0. 

Theorem 2 Let (M, g) be a solution of the Einstein equations coupled to a 
scalar field in four dimensions with symmetry of Bianchi class A. If the Bianchi 
type is I, II or VIq then the variables in the Wainwright-Hsu system converge 
to limits as t — t 00 . If the Bianchi type is VIII the u-limit set of each solution 
is empty. In type VHq both types of behaviour occur. The sets of points which 
occur as uj -limit points is the disk T> in type I, the open subset of T> defined by 
the inequality E + > | in type II, the point off) with (£+, E_) = (—1,0) in type 
VIq and C in LRS type VHq . The limit of f2 is non-zero in types I and II and 
zero in types VIq and VIII. 

Proof The Bianchi type I solutions of the Wainwright-Hsu system with a stiff 
fluid are time independent and are exactly the points of T>. It follows from 
Theorem [1] that any solution of type II converges to a point of T> as r —> 00. It 
is a consequence of the evolution equation for JVi that this point must satisfy 
E + > 0. A more detailed analysis reveals that this inequality is strict. To 
see this note that all terms on the right hand side of the equations for Bianchi 
type II solutions have a common factor N\. Omitting this factor leads to a new 
system whose integral curves are the same as those of the original system in the 
region where Ni ^ 0. The new system is 



Ni 



2(1 - 2E+) 




(30) 




(31) 



(32) 
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This system has a line of stationary solutions defined by the conditions N\ = 
and £+ = i. Linearizing the system about one of these points reveals that it 
has the eigenvalues 0,±3i. Because of the reduction theorem (cf. [TS], section 
5.6) the qualitative behaviour of the solutions can be determined by analysing 
the dynamics on a centre manifold. The system on the centre manifold is such 
that no solution can converge to the stationary point while staying in the region 
Ni > (cf. [16], p. 144, Theorem 5) and this implies the desired result. Any 
point of V with Eo > \ occurs as the w-limit point of a solution of type II since 
the linearization there has a negative eigenvalue. By Theorem[T]any solution of 
type VIo converges to a point of V. As a consequence of the evolution equations 
for N2 and N$ this can only be the point where E + = — 1. In particular 
Q — > and the solution looks asymptotically like a vacuum solution. Any 
locally rotationally symmetric solution of type VIIo has £_ = 0, N2 — N3 and 
a constant value of £+. As t — > 00 the quantity N2 tends to infinity. For type 
VIIo solutions which are not LRS it has already been shown that Y? + + £1 — > 1 
and that £_ and N2 — N3 tend to zero. If the solution has no uj- limit point 
then fl — > as r — > 00 and E + — > 1. For general type VIII it is shown in [T2] 
that fl — > and N + — > 00. 



4 A criterion for geodesic completeness 

The aim of this section is to prove the following criterion for geodesic com- 
pleteness. The general set-up is a spatially homogeneous spacetime defined 
by a one-parameter family g a b(t) of left-invariant metrics on a Lie group G of 
dimension n. 

Theorem 3 Let (M, g) be a locally spatially homogeneous vacuum spacetime 
of dimension n + 1 < 10 whose symmetry is defined by a Lie group G which 
admits no left-invariant metrics of positive scalar curvature. If for some to the 
generalized Kasner exponents satisfy 

9 — n 

supminp a (£) — inf minp a (t) < (33) 

t>to a *>*o a An 

then the spacetime is geodesically complete. 

Proof Since the scalar curvature is non-positive it follows from the Hamiltonian 
constraint that k a \ 3 k ah < (trfc) 2 and hence that J2 a Pa — !• The minimum of 
J2 a Pa subject to the constraint J2 a P a ~ 1 occurs when all p a are equal, as 
follows from the method of Lagrange multipliers. Hence 

k ab k ab > M!. (34) 
n 

The minimum possible value of p a is — as nas been proved in [18j . p. 96. By 
symmetry this is also the minimum value of any other p a . It can be concluded 
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that the maximum eigenvalue A max satisfies 

1 n — 2 

— trfc < A max < trfc. (35) 

n n 

If the value of A max is fixed the inequality (|34l) can be improved. Assume that 
Amax = Ai. Using Lagrange multipliers again it can be seen that the minimum 
of k ab k ab is attained when all A a other than Ai are equal. Hence the lower bound 
for k a bk ah is got by setting 

A 2 = ... = A w = trfc ~ A 1 m " - ( 36 ) 
n—1 

It follows that k ab k ab is bounded below by 
A 



2 (trfc A max ) 2 



max 



n- 1 

= (trk) 2 (np 2 min - 2 Pmin + l)(n - l)" 1 

- (trfc) 2 P(p min ). (37) 

The function P is positive and decreasing on [— 2=2, 11 , Let po be the supre- 
mum in t of p m i n (t). Then by definition p min (t) < po for all t. On the 
other hand the hypotheses of the theorem imply that p m in(i) > Pa — A where 
A = po — inf t>t maL a pa(t). Assume that po — A < 0. (Geodesic completeness 
in the case po > A is treated in Theorem HJ) Now 

d t (trk) = k ab k ab > P(p )(tvk) 2 . (38) 

Integrating this in t gives 

(trfe)(i) > "^y/ (39) 

for a constant C depending on to and trfe(to). Let q a be the projection of the 
tangent vector to a causal geodesic onto the hypersurfaces of constant t. Then 

J t (9a b q a q b ) = 2k ab q a q b < 2(p - A)(trk)g ab q a q b . (40) 
Replacing trfc in this inequality using (1391) and integrating gives 

{g ab q a q b r h >C{l + t)- x (41) 

where A = • The affine parameter length of a causal geodesic up to time r 
is given by 

C (e + g ab q a q b )-Ut (42) 

where e is one for timelike geodesies parametrized by arc length and zero for 
null geodesies. This integral diverges provided A < 1. It remains to show that 
this follows from the inequality A < ^f 1 - Now 
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is equivalent to 

, n — 3 1 ra — 1 A „ ,,,, 

+ Po + A > 0. 44 

n n n 

The expression on the left hand side of this inequality has its minimum at 
Po = —^r- Its value there is 

n — 1 / 9 — n 



4?? 

This completes the proof of the theorem. 



-A . (45) 



For simplicity this theorem has been stated for the vacuum case only. In 
fact the theorem, and its proof, generalize straightforwardly to the case where 
matter is present provided the matter satisfies the dominant and strong energy 
conditions and has reasonable evolution properties. The precise formulation of 
the last property is the matter continuation criterion (MCC) introduced in [18]. 
A variant of Theorem [3] is given by 

Theorem 4 Let (M, g) be a locally spatially homogeneous vacuum spacetime of 
dimension n + 1 whose symmetry is defined by a Lie group G which admits no 
left-invariant metrics of positive scalar curvature. If for some to the generalized 
Kasner exponents satisfy p a ~>—^ then the spacetime is geodesically complete. 

Proof The proof proceeds in a similar way to that of the previous theorem. 
This time the inequality dtittk) > ^ trfc -* implies that 

(trfc)(t)>-^ (46) 

and that 

(WY)" 1 >C(l+t)-\ (47) 

This completes the proof. 

This theorem also generalizes easily to the case with matter satisfying the DEC, 
SEC and MCC. In the case n = 3 the hypothesis on the p a is satisfied identically. 

Next a large class of Lie groups will be exhibited for which the assumptions of 
Theorem[3]are satisfied. They are obtained using products of three-dimensional 
Lie groups with Abelian groups of arbitrary dimensions. 

Theorem 5 Let (M, g) be a spatially homogeneous solution of the vacuum Ein- 
stein equations in n+1 dimensions with 3 < n < 9. Suppose that the group defin- 
ing the homogeneity is of the form G3 x (M x Z2) ra_3 for some three-dimensional 
Lie group G whose Bianchi type is neither IX or VI _i. If the solution is max- 
imally extended to the future then it is geodesically complete. 

Proof Under the hypotheses of this theorem Kaluza-Klein reduction can be 
used to relate the generalized Kasner exponents of the higher dimensional vac- 
uum model to the generalized Kasner exponents of the reduced Einstein-scalar 
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field model in four dimensions. By Theorem Q] and the formulae (fT5T) and (fT¥|) it 
follows that the generalized Kasner exponents of the (n + l)-dimensional space- 
time converge as t — > oo. This conclusion makes use of the fact that, since 
= (^y) 2 , the convergence of ^ is equivalent to that of fi. Hence po = 
and the hypotheses of Theorem [3] are satisfied. 

5 Further examples 

In the previous section geodesic completeness was proved for models based on 
Lie groups belonging to a certain class not admitting metrics of positive scalar 
curvature. To go further it is necessary to enter into the structure of more 
general Lie algebras. (A good general reference for the theory of Lie groups 
and Lie algebras is [22].) For this some definitions are necessary. If L is a Lie 
algebra denote by [L, L] the subalgebra generated by the elements of L of the 
form [x, y] for x and y in L. Define recursively by L^ n+1 ^ — [D- n > , L^] and 
L^> = L. If £(") = {0} for some n then L is called solvable. The radical R of 
a Lie algebra L is an ideal of L which is solvable and which contains all other 
solvable ideals. A Levi subalgebra is a subalgebra which, as a vector space, is 
complementary to the radical. It is semisimple which means that it has no non- 
trivial solvable ideals. A derivation of L is a linear mapping D which satisfies 
D([x,y]) = [Dx,y] + [x, Dy] for all x and y. A particular type of derivations 
are the inner derivations which are of the form [x, ] for some element x of the 
algebra. Let L\ and L 2 be two Lie algebras. Denote by DerL2 the vector space 
of all derivations of L 2 ■ It is closed under commutators and in this way acquires 
the structure of a Lie algebra. Let be a Lie algebra homomorphism from L\ 
to DerL2- A Lie algebra L can be defined as the direct sum of the underlying 
vector spaces of L\ and L 2 with the Lie bracket defined by 

[{xi,yi), (2:2,2/2)] = ([£1,2:2], [2/1,2/2] + 4>[x 1 )(y 2 ) - (j)(x 2 )(yi)). (48) 

The Lie algebra L is called the semidirect sum of L\ and L 2 - The Levi-Malcev 
theorem says that any finite-dimensional Lie algebra can be written as a semidi- 
rect sum of its radical R with a semisimple Lie algebra. 

Consider now the case of four-dimensional Lie algebras. It is known that 
the only semisimple Lie algebras of dimension no greater than four are su(2) 
and sZ(2,R). It follows that the only four-dimensional Lie algebras which are 
not solvable are semidirect sums of the real numbers with su(2) and sl(2,R). 
In fact the semidirect sum of the real numbers with a semisimple Lie algebra is 
isomorphic to a direct sum. To see this note that every derivation of a semisimple 
Lie algebra is inner. Hence the mapping <j) defining the semidirect sum must 
be of the form <p{x) = 4>s(x)[ya, ] for & linear map 4>s from L\ to R and some 
Vo £ £2- It then follows that the linear map ^ from L to itself defined by 

i/j(x,y) = (x,y + (j>s{x)y Q ) (49) 

is an isomorphism from the given Lie algebra to the direct sum. The group 
SU{2) x K admits left-invariant metrics of positive scalar curvature because 
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SU(2) does. On the other hand it is well-known that the simply connected 
Lie group with Lie algebra sZ(2,R) is diffeomorphic to R 3 . Hence the group 
corresponding to the Lie algebra s/(2, R)ffiR is diffeomorphic to ]R 4 . Any simply 
connected solvable Lie group is diffeomorphic to R™ for some n and hence admits 
no metric of positive scalar curvature. Putting these facts together it can be 
concluded that the only connected and simply connected four-dimensional Lie 
group which admits a metric of positive scalar curvature is SU(2) x R. This 
fact was noted in the thesis of Hervik [S], who refers to a paper of Patera et. 
al. [TS]. In the latter a list of Lie algebras of dimensions 4 and 5 is given 
which are not direct sums and they are all stated to be solvable. It may be 
noted in passing that on the basis of this the only types for which recollapse 
can occur in five space dimensions are su(2) © R 2 and su(2) © A 2j \. A general 
criterion for the existence of left- invariant metrics of positive scalar curvature 
in any dimension has been given in [10j . The facts just mentioned suggest the 
following conjecture: 

Conjecture Let G be a connected and simply connected four-dimensional Lie 
group which is not isomorphic to SU (2) x R. Then any spatially homogeneous 
solution of the vacuum Einstein equations with this symmetry group which is 
expanding at some time is future geodesically complete. 

Theorem [5] proves this conjecture for Lie algebras which are the direct sum 
of a three-dimensional Lie algebra and the real numbers and metrics which have 
an additional discrete symmetry. Another case where something is known is 
that where the four-dimensional algebra is a direct sum of two non-Abelian 
subalgebras |Tj. In the notation of [TS] and [5] this is the case ^2,1 © ^b.i- The 
metrics covered by the theorem of [T] are not the most general left-invariant 
metrics. They have a six-dimensional symmetry group. By Theorem 3.1 of [lj 
the quantities P and Q converge to positive limits as t — > 00 and these are 
essentially the generalized Kasner exponents. Hence Theorem @] above applies. 
Similar results are obtained in higher dimensions. These apply for example to 
Lie algebras of the form A 2> i © R n for any n. The limits of the generalized 
Kasner exponents are non-negative. It is once again the case that only metrics 
with additional symmetry are included. 

The Lie algebras which define Bianchi models of class A are those which are 
unimodular, i.e. those for which the trace of the structure constants is zero. 
A four-dimensional Lie algebra which is the direct sum of a three-dimensional 
algebra with the real numbers is unimodular if and only if the three-dimensional 
Lie algebra is of class A. The unimodular four-dimensional Lie algebras which 
are not of this form can be read off from Table 1 in [5]. In the notation used 
there they are A^x, A^, A pq 5 for p + q = — 1, A pq 6 for p + 2q — 0, A^^ 
and A^io- Some explicit solutions are known. In the cases A^ 2 and ^4^5 
explicit diagonal solutions are given in [5] . In both cases the generalized Kasner 
exponents converge to (1, 0, 0, 0) and thus these solutions are future geodesically 
complete by Theorem 0] In the case A^g a solution which is diagonal and self- 
similar is given in [3J. Its generalized Kasner exponents are time-independent. 
They are equal to (§,— \, |, |)- Theorem [4] applies to this solution and it is 
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interesting to note that this is a borderline case for that theorem. In the case 
^4,10 an explicit solution is mentioned in [3] but the authors note that it has a 
higher dimensional symmetry group and can also be considered as corresponding 
to the Lie algebra which is the direct sum of the Bianchi type II Lie algebra 
with the real numbers. It admits a reflection symmetry of the type which makes 
it belong to the class covered by Theorem [5] 

There is a way of fitting many of these Lie algebras in a common geometrical 
framework. In [T7] the notion of local U(l)xU(l) symmetry was introduced and 
it was shown how the Bianchi types I, II, VIo and VIIo can be identified among 
spacetimes with two commuting spacelike local Killing vectors and a global 
topology which is a bundle over the circle with fibre T 2 . The different Bianchi 
types are distinguished by the algebraic structure of a certain 2x2 matrix with 
unit determinant. If both eigenvalues are real with unit modulus and the matrix 
is diagonalizable then type I is obtained. A matrix which is not diagonalizablc 
gives rise to type II. When the eigenvalues are real and distinct type VIo results. 
Finally type VIIo is obtained when the eigenvalues are imaginary. The same 
analysis can be repeated for one dimension higher, although it does not cover 
all relevant algebras. In this case local U(l) x 17(1) x U(l) symmetry is the 
basic input and different cases are classified by the algebraic properties of a 3 x 3 
matrix with unit determinant. If there is a real eigenvalue with unit modulus and 
the corresponding Jordan block is one by one then a case is obtained where the 
Lie algebra is the direct sum of one of the algebras seen in the three-dimensional 
case with the real numbers. The only other case where all eigenvalues have are 
real with unit modulus is that of a 3 x 3 Jordan block. This gives the algebra 
Ai t i. When the eigenvalues are real and all distinct type A^ q 5 results. In the 
case that there is pair of non-real eigenvlaues the Lie algebra which arises is A* q 6 . 
If there is a non-trivial Jordan block corresponding to a double eigenvalue then 
A^ 2 occurs^. This construction does not apply to types A4.8 and A4.10 since 
it can be shown straightforwardly that neither of these has a three-dimensional 
Abelian subalgebra. This point of view shows that in all the cases where the 
Lie algebra has a three-dimensional Abelian subalgebra the group can be made 
to act on the universal cover of a spacetime with a compact Cauchy surface by 
isometries of the pull-back of the spacetime metric. The topology of the Cauchy 
surface is that of a bundle over S 1 with fibre T 3 . 

A convenient feature of vacuum Bianchi models of class A is that initial 
data whose components in a suitable basis are diagonal lead to solutions which 
are also diagonal in that basis. Call a basis of this kind a canonical basis. 
It is moreover the case that for any initial data of Bianchi class A there is a 
canonical basis in which it is diagonal. These facts together mean that in order 
to understand the dynamics of vacuum models of Bianchi class A it is enough 
to understand the dynamics of diagonal models. The fact that models which 
are initially diagonal remain diagonal can be proved using automorphisms of 
the Lie algebras involved. For convenience let the linear transformation which 

1 Remarks related to the material of this paragraph are contained in an unpublished 
manuscript of S. Hervik 
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multiplies by (— l) ai , ai <G {0,1}, be denoted by (010203). In a canonical 
basis the transformations (110), (101) and (Oil) are automorphisms. In other 
words, they leave the structure constants invariant. Data which are diagonal in 
the given basis are invariant under these transformations. Because the structure 
constants are invariant the Einstein evolution equations are also invariant under 
these transformations. It follows from the unique determination of solutions by 
data that invariant data lead to invariant solutions. A solution which is invariant 
under all these transformations is diagonal. It would be nice if this argument 
could be generalized to the case of unimodular Lie algebras in higher dimensions. 
This is no problem if the Lie algebra is a direct sum of a three-dimensional Lie 
algebra with the real numbers. For in that case the automorphisms of the 
three-dimensional case can be trivially extended to get (in a hopefully obvious 
notation) the automorphisms (1100), (1010) and (0110). For the other Lie 
algebras this method has very limited success, at least if the bases in which the 
Lie algebras are presented in [8] are taken as candidates for a canonical basis. 
The only case in which there are enough automorphisms obtained by changing 
the signs of some of the basis vectors to make the argument go through is . 
For that Lie algebra any automorphism must fix en but the signs of the other 
basis vectors can be changed independently. In other words it is possible to use 
the automorphisms (1000), (0100) and (0010). The Wainwright-Hsu system has 
been very useful for the study of Bianchi class A models. Its simplicity relies 
on the diagonalization property and so it may be that nothing of comparable 
power is available for general unimodular Lie algebras in any higher dimension. 

The theme of the present paper is the late-time behaviour of cosmological 
models but the behaviour near the initial singularity is also of great interest. 
For the latter subject the issue of diagonalization plays a very important role. 
In spacetime dimension four it is believed that generic spatially homogeneous 
solutions of the Einstein vacuum equations show oscillatory (Mixmaster) be- 
haviour near the initial singularity and this has been proved rigorously in some 
cases. (For Bianchi type IX see [T5].) In 4 + 1 dimensions it seems that this is 
still true but that restricting to the special case of diagonal models removes the 
generic oscillations [B|. In particular, generic oscillations can occur for the group 
SU{2) x K due to the influence of the non-diagonal degrees of freedom. The 
situation is similar in dimensions up to 9 + 1 but beyond that the oscillations 
are no longer generic. For some rigorous results on this see [4]. 

6 Conclusions and outlook 

In the preceding sections it has been proved that spatially homogeneous solu- 
tions of the Einstein equations in various dimensions with symmetries defined 
by a large class of Lie groups are geodesically complete in the future. The case 
considered in most detail is that of spacetime dimension five. Even there it 
is to be hoped that something much more general can be proved. A sugges- 
tion for a generalization of this kind is formulated as a conjecture in Section 
5. Making a corresponding conjecture in general dimensions is dependent on 
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a sufficently good understanding of the structure of general Lie algebras. A 
guide in formulating the conjecture was to restrict to those Lie groups which do 
not admit a metric of positive scalar curvature. Another direction of possible 
generalization is to consider spatially homogeneous spacetimes where the Lie 
group defining the symmetry is necessarily of a dimension higher than the space 
dimension. These generalize the Kantowski-Sachs models in spacetime dimen- 
sion four. In five spacetime dimensions the list of possibilities is as follows [8]: 
the four-sphere 5 4 with symmetry group 5*0(4), the complex projective plane 
CP 2 with an eight-dimensional symmetry group, S 2 x S 2 with a six-dimensional 
symmetry group, S 2 x R 2 and S 2 x H 2 with five-dimensional symmetry groups. 
All of these examples admit an invariant metric with positive scalar curvature 
and so are not promising candidates for allowing spacetimes which are future 
geodesically complete. The theorem of Berard-Bergery [2] mentioned at the end 
of Section 2 extends to homogeneous spaces which are more general than Lie 
groups and may be helpful in investigating generalizations of Kantowski-Sachs 
models in higher dimensions. 

In the previous sections the cosmological constant A has always been as- 
sumed to vanish. It is likely that the theorem of Wald [24] on the late-time 
asymptotics of spatially homogeneous cosmological models with A > extends 
to arbitrary dimensions with the essential assumption (in the vacuum case) be- 
ing the absence of metrics of positive scalar curvature. Combining this with the 
ideas in [T3] should give geodesic completeness. These ideas have, however, not 
yet been worked out in detail. 

There are very few results on future geodesic completeness for inhomoge- 
neous spacetimes which do not rely on choosing initial data which are close to 
those of a known solution (small data results) . This applies even in the a priori 
simpler case A > 0. In fact the only result of this kind for vacuum spacetimes 
is that of Ringstrom [21]. It could be that obtaining a better understanding 
of geodesic completeness for homogeneous models in higher dimensions would 
provide new ideas for understanding the inhomogeneous case in spacetime di- 
mension four. 

It is interesting to ask what happens when there is a metric of positive scalar 
curvature. Is there always recollapse in that case? It seems that very little is 
known about this in spacetime dimensions greater than four. There are also 
few results available for inhomogeneous solutions in four dimensions. Perhaps a 
better understanding of the homogeneous case in higher dimensions could help 
to change this. In particular, higher dimensions could be a good place to look 
for a counterexample. 

The long-time behaviour of solutions of the vacuum Einstein equations in 
higher dimensions gives rise to a rich variety of problems concerning the dy- 
namics of the gravitational field. In this paper we have obtained theorems on 
some of the simpler of these, while pointing out promising directions for future 
research. We have also emphasized that one motivation for studying these ques- 
tions is the insights they may bring for the Einstein-matter equations in four 
dimensions. 
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